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Abstract. Motivated by recent experimental results and ongoing measurements, we review the chiral per- 
turbation theory prediction for Kl — * it^e decays. Special emphasis is given to the stability of the 
inner bremsstrahlung-dominated relative branching ratio versus the K e % form factors, and on the sep- 
aration of the structure-dependent amplitude in differential distributions over the phase space. For the 
structure-dependent terms, an assessment of the order p 6 corrections is given, in particular, a full next- 
to-leading order calculation of the axial component is performed. The experimental analysis of the photon 
energy spectrum is discussed, and other potentially useful distributions are introduced. 

PACS. 13.20.Eb Radiative semileptonic decays of K mesons - 11.30.Rd Chiral symmetries - 12.39.Fe 
Chiral Lagrangians 
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1 Introduction 

The amplitude for the semileptonic radiative decays Kl — > 
7r =F / ± ^7 [-Kj3 7 ] , with I = e,fx, can be divided into two com- 
ponents: the inner bremsstrahlung (IB) that accounts for 
photon radiation from the external charged particles and 
which is determined by the non-radiative process Kl — * 
tt^^vi [if/3]; and the structure-dependent (SD) ampli- 
tude, also called "direct emission", that describes pho- 
ton radiation from intermediate hadronic states and rep- 
resents genuinely new information with respect to the IB 
one. 

Low's theorem pQ, applied to if/3 7 , states that the 
leading contributions in the expansion of the amplitude 



in powers of the photon four-momentum q, namely, the 
orders g _1 and q° , are completely determined in a model- 
independent way by the IB via the if/3 form factors and 
their first order derivatives. The SD amplitude is then de- 
fined by the terms of order q and higher. In |21EI], the 
procedure of was followed to derive the g _1 and q° 
terms of the IB amplitudes for if/3 7 ; moreover, a qualita- 
tive, model-dependent, assessment of the SD amplitudes 
was performed using vector meson dominance. In jjj, the 
radiative decay modes for both charged and neutral kaons 
were calculated, taking into account IB terms only. Orig- 
inally, the main interest was the precision test of soft- 
photon theorems, allowed by the dominance of IB and the 
fact that, for if/3 7 , the non-radiative amplitude could in 
principle be studied extensively and with high accuracy. 

Later, if/3 7 decay amplitudes (including charged kaon 
ones) were calculated at order p* in chiral perturbation 
theory (ChPT) in 0, and branching ratios were evalu- 
ated for I = e, 11 in a feasibility study for DAFNE [Hj. An 
error analysis and a dedicated study of decay distributions 
was postponed to a later stage when precise data would 
become available. It is one of the aims of the present work 
to provide such an analysis. 

ChPT allows for a systematic expansion of transition 
amplitudes for low momenta of the pseudoscalar mesons 
[7||H]. The lowest order amplitude is only of the IB-type 
with constant if/3 form factors, and is independent of free 
parameters. In addition to providing momentum depen- 
dence of the Kis form factors in the bremsstrahlung, the 
0(p 4 ) terms predict the existence of no n- vanishing SD am- 
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plitudes (vector and axial- vector) , unambiguously calcula- 
ble in terms of loop diagrams, low-energy constants of the 
strong Lagrangian £4 8 , and the chiral anomaly [9II10| . 
While the anomaly does not require new physical param- 
eters, the low-energy constants are numerically already 
well-determined from other, independent, meson pro- 
cesses. Consequently, the experimental verification of the 
SD amplitude currently represents a significant test of 
ChPT and, ultimately, of QCD. Of course, since the ex- 
pansion of the SD amplitudes starts at 0(p 4 ), one may 
inquire about the role of higher order corrections, a point 
that will be addressed in the sequel. 

In practice, this experimental analysis is complicated 
by the fact that the radiative K13 branching ratio is largely 
dominated by the IB, while the SD contribution via IB- 
SD interference is expected to be an effect at the percent- 
level (the pure SD rate is negligibly small). On the other 
hand, the characteristic q^ 1 behavior of the IB by far 
dominates the lower (and intermediate) photon energy 
range 0], while in the upper range where the SD effects 
become more significant, the number of events is severely 
reduced. In addition, precise knowledge of the K13 form 
factors is required for a reliable fit to the photon energy 
distribution, in order to improve the sensitivity to signals 
of SD contributions through deviations from the pure IB. 
All that calls for high precision measurements of both 
and K13. 



The first measurement of the decay Kl 



with significant statistics was performed by the NA31 Col- 
laboration in a pioneering experiment, which proved 
the possibility of precision measurements of this process. 
Their result for the decay rate relative to K e % decays and 
for the branching ratio agreed with the theoretical pre- 
dictions of @] and [5], respectively. A few years later, the 
KTeV Collaboration ^2] determined the relative branch- 
ing ratio, together with the photon spectrum, at percent- 
level sensitivity. They found a result which is "significantly 
lower than all published theoretical predictions" . More- 
over, from the measured photon energy distribution, two 
particular combinations of the SD amplitudes were ob- 
tained, within rather large uncertainties. Quite recently, 
new experimental results on K e 3j with percent-level ac- 
curacy have been presented by the NA48 ^21 an d by the 
KTeV ^2] Collaborations, respectively. 

The experimental situation has therefore become quite 
promising and clearly justifies renewed interest in radia- 
tive K e 3 decays. In this regard, particularly relevant top- 
ics are the stability of IB with respect to the K e 3 form 
factor parameterizations, the separation itself of the ra- 
diative amplitude into IB and SD contributions and the 
re-visitation the ChPT calculation of the SD component, 
including in particular next-to-leading 0(p 6 ) corrections. 
Moreover, this represents an opportunity to discuss, be- 
sides the photon spectrum, also other differential distribu- 
tions that may help in the study of the SD terms in high 
statistics experiments. 

In the sequel, we limit our consideration to the K e ^ 
transition, since current experimental data with appro- 
priate statistics refer to this mode only. Specifically, in 



Sects. and El we review the experimental situation and 
the experimental observables, the kinematics, and the am- 
plitude definitions with particular emphasis on the sepa- 
ration into IB and SD contributions. In Sect. 01 we present 
the ChPT results for the SD amplitudes, notably the 
0(p 6 ) corrections for the axial amplitudes. In Sect. we 
numerically discuss the relative radiative branching ratio, 
while Sects. El and [7| are devoted to numerical estimates of 
the SD amplitudes, the photon energy distribution and the 
comparison with experimental results. Also, other kinds 
of differential distributions are discussed there. Finally, 
Sect. El contains a summary of the results, while details of 
the calculations are collected in the appendices. 



2 Experimental status and observables 

As already mentioned, we concentrate on decays 
where high statistics experimental data on the branching 
ratio and photon energy spectrum have become available. 
For a presentation of the experimental status in the other 
channels, we refer the reader to the PDG listing |15| . 

Measuring the decay rate relative to K e 3 is much safer 
than absolute measurements, as the former is free from 
uncertainties related to experimental normalizations, cal- 
ibrations, and machine luminosity. Basically, an inclusive 
K e 3 sample of events is collected, all characterized by one 
lepton and one pion of opposite charges emitted from a 
common vertex, without any restriction on the number of 
emitted photons. A radiative K e 3 7 subsample is extracted 
by imposing additional criteria dictated by the apparatus 
and the experimental conditions, in particular by the re- 
quest of having at least one hard photon in each of those 
events. 

To achieve optimal identification of the candidate K e 3y 
events, kincmatical cuts are applied to the radiative sam- 
ple. In particular, thresholds in the photon energy and in 
the photon-electron opening angle are usually imposed, 
see, e.g., 0]. Then, the experimental results generally con- 
cern the relative branching ratio 



r(K, 



J?* \ 77'CUt /)* \ /lCUt 



r(K e3 ) 



_ ^ (2.1) 

where E* and #* 7 indicate the photon energy and the 
photon-electron opening angle in the kaon rest frame, re- 
spectively. 

From the above, the measured value of R is deter- 
mined by the ratio of the number of events in the K e3l 
and K e 3 samples, each divided by the respective exper- 
imental acceptances. The available experimental results 
are displayed in Tabled 

One may notice that the KTeV 04 result for R [HI with 
angular cuts is based on a much smaller number of K e 3j 
events with respect to their previous determination [12j . 
yet the achieved uncertainty based on their more recent 
data is comparable owing to a much reduced systematic 
uncertainty. 

As regards the photon energy distribution, one can in- 
vestigate the spectrum with free normalization because 
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Table 1. Experimental values of R for the transition K e z-y. 
The first error is statistical, the second one systematic. 



Ref. E% 



events 



R x 10 2 





30 MeV 


20° 


18977 


0.964 ± 0.008 t ° °ool 


M 


30 MeV 


20° 


4309 


0.916 ±0.017 


m 


10 MeV 


0° 


14221 


4.942 ±0.042 ±0.046 


m 


30 MeV 


20° 


15463 


0.908 ± 0.008 t 0^12 


cu 


30 MeV 


20° 


1384 


0.934 ± 0.036 ± ° ° ll 


M 


15 MeV 


0° 


10 


3.3 ±2.0 



the essential features lie exclusively in the shape. Indeed, 
structure-dependent emission should manifest itself in the 
harder portion of the photon energy spectrum, via a modi- 
fication of the pure IB spectrum which is controlled by the 
K e 3 form factors. An attempt along these lines was made 
by the KTeV 01 experiment |12| . which used a simpli- 
fied decomposition of the structure-dependent amplitude, 
and their analysis will be commented upon in the sequel. 
This is the only experimental information on structure- 
dependent emission currently available, as neither NA48 
nor KTeV have so far presented an analysis of this point 
based on the more recent data. 

As far as the perspectives of measurements are 

concerned that may complement the currently available 
data for K,,^, new results on these transitions are ex- 
pected from the NA48 experiment |17| . In the more remote 
future, substantially increased statistics for decays 
should be expected in connection with the construction of 
higher intensity kaon beams, with accumulated samples of 
10 6 (or more) candidate events [T5|. 



3 The decay amplitude 

In the following, we consider the decay 

K\p)^ir-(p')e + {p e ) V M 1 {q) [X° 37 ] (3.1) 

and its charge conjugate mode. We disregard CP-violating 
contributions, and study the emission of a real photon 
(q 2 = 0). 

3.1 The matrix element 

The transition matrix element has the form 



T(K°) = % e y u * s e*(g) 



V2 



x u(p v ) 7 W (1 - 75) v(p e ) 

V>(P„) I" (1 ~ 75) (m e - i>e- i) 7m v (Pe) 



W, e 





a) Njc" b) 
Fig. 1. Diagrams describing K° 3l decay 

The relevant diagrams are displayed in Fig. The first 
term of l|3.2[) corresponds to diagram a), which includes 
bremsstrahlung off the charged pion, while the second one 
corresponds to the radiation off the positron, represented 
by the diagram b) . We have introduced the hadronic ten- 
sors Vpv and 

I^=i J <£xe^ (n-(p')\TV™(x)I^(0)\K°(p)) ; 
I = V,A, 

(3.3) 

whereas is the K13 matrix element 

Ffi = (■K-(p')\V^(0)\K°(p)) , (3.4) 



with 



^ ad = 5Tp« , = s "^75« 



V° m = (2uu-dd- ss)/3 . 
The tensors and satisfy the Ward identities 
<fV* v = F v , 



(3.5) 



<fA„ v = 



(3.6) 



(3.2) 



which imply gauge invariance of the total amplitude (|3.2I) , 
<f M p = . (3.7) 



3.2 Inner bremsstrahlung, structure-dependent 
terms and all that 

Low's theorem is employed in [2] to obtain the IB am- 
plitude for Ki3 7 decays, written entirely in terms of the 
K13 form factors and their derivatives. In this subsection 
we present an alternative way of separating the amplitude 
into an IB and a SD part that directly starts from the 
following two requirements: 

1. In order to describe two different physical mechanisms, 
the IB and SD amplitudes must be separately gauge 
invariant. 

2. The SD amplitude contains terms of order q and 
higher. 

The second condition does not prevent the IB amplitude 
from containing terms of order q and higher. Splitting the 
amplitude under this less restrictive condition allows one 
to put more terms into the IB part, still using only the 
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non-radiative matrix element in this part of the amplitude. 
This has the advantage that one can make more precise 
predictions for the decay process, as we will see below. 

The splitting of the transition amplitude T into an IB 
and a SD part requires a corresponding splitting of the 
hadronic tensor I^ v . Consider first the axial correlator. 
There are no contributions where the photon is emitted 
from the pion line. Therefore, is considered to be a 
purely SD contribution. It can be written in the form 5 

A - A SD 

= % e^pa {A 1 p'Pq* + A 2 q(>W a ) (3.8) 
+ i e M Apa p' x q p W° (A 3 W v + Ai v ' v ) , 

where = (p — p' — (We use the convention eom = 
+ 1.) The amplitude is manifestly of order q and higher, 
because the Ai are non-singular at zero photon energy. 
The Lorentz invariant components Ai depend on three 
independent scalar variables that can be built from p, p' , 
and q - we come back on this in the following paragraph. 
The decomposition of the vector correlator reads 



V. — V m + V au 
where the IB piece is chosen such that 

<fV™ = > ( 3 - 10 ) 

as a result of which we have 

g*V* D = . (3.11) 

The structure-dependent part of the decay amplitude T 
in (|3.2() is defined to be 



^SD 



(3.9) 



Gi 
V2 



(3.12) 



x u{p v )Y (1 -7s)«(pj) , 



whereas the bremsstrahlung part is T IB = T — T SD . 

It remains to explicitly construct the decomposition 
<|3.9[l . In order not to interrupt the argument, we refer the 
interested reader to Appendix [E| and simply display here 
the result, 

V™ = ^(2 P MW 2 )-W u f 2 (W 2 )) 
p'q 

+ ^(2( P -q)„Af + -W v Af 2 ) (3.13) 

+ g MV (2Af+-f 2 (t)) , 

A fa = fi(t)-fi(W 2 ) , i = +,2 , (3.14) 

where /+, f 2 are the form factors Ij3.4|l 

F^ = 2p^ + (t) + (p'-p)J 2 (t) , t=[p-p'f . (3.15) 

We use the form factors / + , f 2 instead of the usual / + , 
/-=/+ — f 2 ones for easier comparison with the work 

of El- 



The IB part derived in differs from the one used here 
through terms of order q. It can be obtained from 
by subtracting all terms of order q and higher from the 
latter, and merging them into the SD part of the ampli- 
tude. Because the terms to be subtracted can be expressed 
through the form factors /+, f 2 and their derivatives, we 
believe that it does not make much sense to perform this 
purification of the IB part, and we will mostly stick with 
the convention l|3.13|l . While comparing with the KTeV 
result ^2jj w e will have the occasion to compare Ij3.13|l 
with the conventional decompositions [3] in more detail in 
Sect. IO 1 

Let us shortly discuss the salient features of the IB 
term (|3.13ll . First, it satisfies the Ward identity (|3.10(l . 
Second, it contains all infrared singular pieces propor- 
tional to 1 1 p'q. With this we mean the following. The 
residue of the singularity is a non-trivial function of the 
momenta p, p', q. The decomposition (|3.13(l takes into 
account all singularities at p'q — 0, in contrast to the 
standard treatment |3], which considers e.g. a term like 
(qW) 2 /p'q to be of order q, to be relegated to the SD part 
of the amplitude. 

It is useful to decompose also the SD part of the vector 
amplitude into a set of gauge invariant tensors. In the 
following, we often use the basis proposed in |20|. 



V, 



Vi^qu ~ p'qg^u) + V 2 (W^q u 
+ V 3 (qWp^W u -pqW^W») 
+ V i (qWp\ l p' u -p'qW ll p' v ) . 



qWg^ 



(3.16) 



The Lorentz invariant amplitudes Vi again depend on the 
3 scalars that can be formed from p, p' , and q. 



3.3 Kinematics 

It remains to shortly recall the kinematics of this decay, 
and we begin with the Lorentz invariant amplitudes Aj, V^. 
As already mentioned, these are functions of three scalar 
variables that we often take to be 

a = {q+p'f , t = (p-p'f , u= {p-qf . (3.17) 

These variables are useful in the discussion of the ana- 
lytic properties of Vi, Ai. In l|3.3l) . the variables s, t, u can 
assume any value, whereas the physical region in K e3l 
decays can be represented as follows. For fixed W 2 , the 
variables s, t, and u vary in 



W 2 < t < (M K - M^) 2 , 



s_ < s < s 



+ • 



S± =Ml--(t + Ml - M 2 K )(t - W 2 ) 

± ^ 2 (t,M 2 K ,M^ 2 (t,0,W 2 ) , 
s + t + u = M 2 K + M 2 + W 2 , 



(3.18) 



Our separation into IB and SD contributions is very close in 
spirit to the notion of generalized bremsstrahlung as developed 

in nnj. 
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where 



4 Analytical results from ChPT 



X(x, y, z) — x 2 + y 2 + z 2 — 2(xy + xz + yz) . (3.19) 

Varying the invariant mass squared W 2 of the lepton pair 
in the interval 



m 2 < W 2 < (M K - M T ) 2 



(3.20) 



generates the region covered by s,t,u in decays. In 
Sect. 14.21 where the analytic properties of the amplitudes 
Ai, Vi are discussed, we display the region l|3.18|l in the 
Mandelstam plane. 

Instead of s,t,u, we also use 



pq/M K = El , PP '/M K = E* 



W 



(Pi +Pv) 2 , 



(3.21) 

where E* E* are the photon and the pion energy in the 
kaon rest frame. This set is useful when discussing partial 
decay widths. 

In the case of four body decays we have five indepen- 
dent variables, thus two more variables are needed to de- 
scribe fully the kinematics of Ki 3l decays. We choose 



PPe /M K = E* , x = p e q/M 2 K , 



(3.22) 



where E* is the positron energy in the kaon rest frame. 
The dimensionless variable x is related to the angle 6* 
between the photon and the positron: 



xM'k = EZ [ El - \/E. 



mi cos (9 



(3.23) 



The total decay rate is given by 

r(K° TT-e+v-f) = (3.24) 
2MA)» J dvrsto p',Pe,Pv,q) J2 ^ > 

spins 

where T is the amplitude in (|3.2|l . and we denote the 
Lorentz invariant phase space element for the Ki 3j process 
by dhipsip', P'tPetPu, q) 2 The square of the matrix ele- 
ment <|3.2[l . summed over photon and lepton polarizations, 
is a bilinear form in the invariant amplitudes Vi,Ai,f+ 
and /2. Performing the traces over the spins, we work 
with massless spinors, as a result of which the form factors 
A3, V3 and /2 drop out in the final expressions. [The elec- 
tron mass cannot be set to zero everywhere, because the 
IB part of the transition amplitude contains mass singu- 
larities, generated by the diagram Fig.^.] In AppendixlBl 
we display the explicit result for X] spins 1^1 > m particular 
the T-odd terms that are generated by the imaginary parts 
of the structure functions, and comment on the relation 
to the width of Kl- 

2 For the decay of a particle of momentum p into n particles 
of momenta pi, . . . ,p„, one has 

d 3 p k 



d L iPs(p;pi, ■ ■ ■ ,Pn) = S^ip-^pij Yl 



While Low's theorem furnishes a recipe to evaluate the 
terms of order 1 /q and q° of an amplitude associated with 
a general radiative process, it does not give any insight 
into the terms of order q and higher, that is the SD part. 
A convenient tool to derive expressions for the SD ampli- 
tude is ChPT. For the axial part, ChPT directly gener- 
ates the corresponding amplitude in a series expansion in 
the momenta, the leading contribution is generated by the 
Wess-Zumino-Witten (WZW) term [101 • ^ s f° r the vector 
amplitude, the chiral expansion contains both IB and SD 
terms, hence if one simultaneously evaluates the K13 ma- 
trix element, the decomposition (|3.13|l leads to the chiral 
expansion of the SD term. 



4.1 ChPT results at order p 4 

In Uj|, the chiral expansion was carried out up to 0(p 4 ) for 
the neutral and for the charged decay modes. [A tree-level 
calculation up to this order without the loop contributions 
was performed in [21] .] We do not describe the calculation 
here and refer the interested reader to the original article. 
The result for the SD terms is as follows. For the axial 
amplitude, one has 



Ao 



1 



8tt 2 F 2 ' 



A l =A 3 = Ai = p(p 4 )} . (4.1) 



F is the pion decay constant in the chiral limit. 3 We dis- 
play the result for the vector amplitude V™* in terms of 
the Lorentz invariant form factors Vi, 



V x = V2I2 
qW 



V 2 = ^(h-p'qL 2 + V2A.f, 



(4.2) 



V 4 = 



[0(P 4 )) 



The integrals Ii , f£ are defined as follows. In |S] , the one- 
loop expression for V^ v in the charged decay mode 
is defined in terms of integrals Ii , , explicitly displayed 
there. The quantities Ii, are obtained from Ii, f% by 

1. replacing the arguments (p,p') by — {p',p) in the Ii ; 

2. inserting the appropriate coefficients c\ for Kf 3 ^ listed 
in Table 10 of that reference. 

Note in addition that A/ + in l|4.2|l refers to the chiral 
one-loop representation of this quantity. 

It turns out that the form factors Vi are nearly con- 
stant over the physical phase space. This is due to the 
fact that in the vicinity of the physical phase space, there 
are no singularities at this order in the chiral expansion. 

3 Usually, the meson decay constant in the SU(3) chiral limit 
is denoted by Fo- We refrain from following this convention in 
order to slightly ease the notation. 
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This fact allows one to derive expressions for the SD parts 
that are considerably simpler than the full ones, yet still 
precise enough for our purpose. In a first step, we expand 
these amplitudes in powers of the photon momentum q 
and keep only the leading order term. This amounts to 
setting s = M%,u = M K , as a result of which the Vi 
become functions of t alone. It turns out that all loop 
integrals can be expressed in terms of the standard one- 
loop integral J(t). The resulting expressions are displayed 
in Appendix El An even more drastic simplification re- 
sults when one furthermore sets t = in the simplified 
formulae. The result reads 



(1 



F 2 J 32tt 2 F 2 
+ (l + x - x 2 + x 3 ) 

- (127 - 93x + 2lx 2 



(53 - 25x + 2x 2 



3 

log a; 



2(1 - x) 



V2 = j2 ( l * 



L w ) 

(l + x)(l-x)- 2 



Q4:ir 2 F 2 



logy 
2(1 -x) 



2x log x 
1-x 



0(q,t) 



(1 



x 



32ir 2 F 2 



166 .„ , . , .x 2 

— (9-4x) + (77-x) — 



V 3 



-x(3 + 2x) 
-0(q,t) , 

32ir 2 F 2 M K 



logx 



l-x 



9(12 -x) (4 -x) 



2 lo S2/ 



2611 . , 4 , 

13a; 34- 5x - -x 3 

3 v ; 3 



-I- x(2 + 3x + x z ) 
+ 0(q,t) . 



2 \0gX 



27(7 - x) (4 -x) 



2 lo S2/ 



(4.3) 



M 2 JM 2 K , y = M 2 /M 2 K . 
derstood that M 2 is related to M K , M 2 through the Gell- 
Mann-Okubo relation. 



4.2 ChPT results at order p 6 

There are two main reasons to consider 0(p e ) corrections 
to the structure-dependent terms as described in the pre- 
vious subsection: 

1. As the structure-dependent terms vanish at tree level 
in the chiral expansion, the one-loop or 0(p A ) predic- 
tions are only the leading order results for these am- 
plitudes. Subleading corrections are often sizeable in 
chiral SU(3), therefore it is mandatory to investigate 
0(p 6 ) terms in order to be sure to control the size of 
the structure-dependent terms. Furthermore, several 
of the structure functions vanish at leading (one-loop) 
order (V4, A\, A3, A4) or nearly so in the sense that 
they do not allow for natural-size counterterms (V3), so 
the size of corrections to these is completely unknown 
from the one-loop results. 




s = s = 9M% 

t = (M K +M n ) 2 




t = 



u = (M K +M n ) 

Fig. 2. Cuts in the complex plane for the vector (left panel) 
and axial (right panel) amplitudes for fixed W 2 = m 2 .. The ar- 
rows indicate the part of the plane where the decay amplitude 
is complex, starting from the lowest cuts indicated by the full 
lines. The limits of physical phase space are also shown (thick 
solid line). 



2. At 0(p 4 ), all structure functions are real in the phys- 
ical region, and the cuts in these functions lie far out- 
side the kinematically allowed range. This is the reason 
why they are so smooth and can be approximated to 
such high accuracy by simple polynomials. However, 
this changes at 0(p 6 ), as will be seen below. 

For the following discussion, we again use the Mandelstam 
variables s, t, u. The lowest-lying cuts for the structure 
functions in terms of these three variables are as follows: 

1. For the weak vector current, they start at s t hr = 9M 2 ., 
tthr = (M K + AU) 2 , u thr = (M K + 2M n ) 2 , respectively. 
Only the t-channel cut exists at 0(p 4 ) as the other two 
require three-particle intermediate states and therefore 
occur only at two-loop order. 

2. For the weak axial current, cuts start at Sthr = 4M 2 , 
Uhv = {M K + M n f, uthr = (M K +M,) 2 , respectively. 
All these occur at one-loop order, but are suppressed 
to 0(p e ) as they require an anomalous vertex. 

These cuts are displayed graphically in Fig. [51 where we 
have drawn them in the Mandelstam plane together with 
the allowed phase space for fixed (and minimal) W 2 = 
m 2 , which corresponds to the maximal range in s, t, u. 
While the t- and it-channel cuts lie far outside the physical 
region, we note that the s-channel cuts overlap with it 
(precisely: for W 2 < (M K -3M^) 2 in the vector and W 2 < 
(Mk — 2Mn) 2 in the axial case), such that at least some 
of the structure functions become complex at 0(p 6 ). 

The diagrams with cuts responsible for imaginary 
parts in the physical region are displayed in Fig. to- 
gether with one typical diagram for the i-channel cut ap- 
pearing at 0(p 4 ). Due to the smallness of phase space for 
the three-pion intermediate state, we expect the effect of 
the cut in the vector structure functions to be tiny. 

For the above considerations, we have regarded W 2 as 
a fixed "mass squared" of the lepton-neutrino pair, which 
is of course not true. There are additional cuts in W 2 , the 
lowest one starting at W 2 = (Mk + M T ) 2 present at one 
loop, but sufficiently far outside the physical region, plus 
a pole at W 2 = M K that however only appears in the 
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Fig. 3. Feynman diagrams with cuts for a) s > 4M^, b) 
s > 9M%, c) t > (M K + A'U) 2 , respectively. The filled ver- 
tex denotes a contribution from the anomalous Lagrangian at 
order p 4 . The first two diagrams provide the structure func- 
tions with imaginary parts in the physical region, while the 
last one generates the only cut present in the amplitudes at 
0( P 4 ). 



structure function A3. The pole at s = M% defines the 
bremsstrahlung part and is not present in the structure- 
dependent terms as defined by our convention. 



slowly above threshold: as the leading and next-to-leading 
order amplitudes in the chiral expansion are real, it is 
certainly negligible in the squared matrix element at our 
accuracy. 

We conclude that even for the axial structure func- 
tions, the impact of the various cuts on the behavior in 
the physical region is rather weak. 

4.2.2 Order p 6 corrections to the vector amplitudes 

A complete evaluation of the vector amplitudes at order 
p 6 requires a full two-loop calculation, which is beyond 
the scope of this article. A less ambitious work consists in 
the determination of the leading chiral logarithms at two- 
loop order [221 • ^ s a m ' s t step in this direction, we have 
explicitly calculated the contributions of the form Li x Lj 
at order p 6 . We find that these can all be absorbed into a 
renormalization of the couplings at order p , according to 
F 2 — > F-xFk- This is strictly analogous to the observation 
that the dependence on L4, L5 for the axial terms can be 
absorbed into such a renormalization. We will make use of 
this fact in Sect. El where we provide numerical values for 
the SD terms. In addition, we will give a rough estimate 
of the contributions at order p e , and defer an evaluation 
of the leading logarithms to a later publication [2*5| . 



4.2.1 Complete order p 6 corrections to the axial amplitudes 

We have calculated the complete 0(p 6 ) corrections to the 
axial structure functions Ai, A2, and A4. A3 is always sup- 
pressed by a factor of mf and is therefore disregarded in 
the context of the electron channel. We find the following 
structures: 



A 2 
A, 



1 



Att 2 F v F k 
1 

8n 2 F^F K 



{S 1 (s)+T l (t)+U 1 (u) + X 1 } , (4.4) 



{l + S 2 (s)+T 2 {t) + U 2 {u)+X 2 } 



FttFj 



K 



(4.5) 
(4.6) 



The explicit forms for the various loop functions as well 
as the combinations of low-energy constants entering the 
expressions I|4.4|) - (I4.6|) can be found in Appendix iDl We 
remark that it is strictly necessary to differentiate between 
F„ and Fx at this order only for A 2 . The normalization 
was chosen this way such that any dependence on the low- 
energy constants L4, L5 vanishes in the final result. 

We furthermore note that only A\ has a contribution 
from the s-channel cut and therefore becomes complex in 
parts of the physical region at this order, while S 2 (s) is a 
simple polynomial. It turns out, though, that also in S\(s) 
the standard two-point loop function J ffW (s) from the two- 
pion intermediate state comes with a prefactor of (s — 
AM 2 ) such that the cusp in the real part is smoothed out. 
This is due to the fact that the 7T7t — > 717 rescattering has 
to be a P-wave and is therefore suppressed at threshold. 
For the same reason, the imaginary part also rises very 



5 The ratio R 

A particularly useful quantity to consider for the exam- 
ination of _Ke3 7 decays is the ratio R defined in 1)2.1(1 . 
rather than the absolute width for the K e ^ channel or 
the branching ratio thereof. This is desirable both from 
the experimental and the theoretical point of view for 
rather similar reasons: both experimentally and theoreti- 
cally, certain normalization factors cancel in the ratio (to a 
large extent), and hence the uncertainties ensuing thereof 
are avoided. To present the situation in a more transpar- 
ent way, we shall initially neglect all possible complications 
ensuing from radiative corrections or isospin breaking, and 
shall discuss these afterwards in Sect. 15.41 In order to re- 
mind the reader of this simplification, we shall denote R 
in the absence of real and virtual photon corrections by TZ 
in the following, 

a^TZ = [a^R]^ . (5.1) 

We may decompose TZ according to 1Z = 1Z IB + 1Z SD in 
the following sense: 7Z lB is understood to be 1Z in the limit 
where all structure-dependent terms are omitted, while we 
may then define 7Z SB =7Z — 7Z lB , such that 1Z SI} contains 
in particular also interference terms of bremsstrahlung and 
structure-dependent terms. 

We begin by deriving a simple expression for 1Z. Start- 
ing from (|3.24l) . we may define a quantity SM by 



r{K e 



1 



3~, ) 



2M k {2k) 1 



J d L1PS \T\ 



spins 



8aM 5 K G 2 F \V us \ 2 f s 
(2nf /+(0) 



(5.2) 



*LIPS 



SM 
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Table 2. Coefficients for the K e s phase space integral. 



Table 3. Coefficients for the K e 3^ phase space integral for 
£^ ut = 30 MeV, 0™* = 20°. The errors for &f D are p 6 estimates. 



ao 


ffli 


ai 


03 


«4 














0.09390 


0.3245 


0.4485 


3.092 


6.073 


°o 


&! B 


62 


bW 


64 


65 












1.509 


5.23 


6.92 


14.71 


47.6 


92.3 



The phase space integral J^lips SM such defined is di- 
mensionless and free of (electroweak) coupling constants. 
For the bremsstrahlung part, also /+(0) 2 factors out nat- 
urally, such that all form factors appearing in SM are the 
normalized form factors f+{t) = /+(£)//+ (0). The non- 
radiative width can be written as 



r(K e3 ) = / dydzp{y,z) 

p&.«) - M \^ t " |2 / + (o) 2 A( g ,z)/ + (t) 2 

where y = 2pp e /M 2 { , z = 2pp' jM\, and 

A(y, z) = i(z + y - 1)(1 - y) + r e (4y + 3z - 3) 
- 4rv + r e (rv - r e ) , 



(5.3) 



(5.4) 



with r e = ml/M^, rv = M 2 /M|-. We therefore find for 
72. the following simple expression, 



K = 



80 



/ d L ip S SM 



7T 4 JdydzA(y,z)f + (ty 



(5.5) 



in which all factors Gf, Kj S , /+(0), and Mk have can- 
celed. [For the relation between K e3l and K® 3j decays, 
see Appendix IBl] 



5.1 Phase space integrals 

Assuming 



t t 2 

f+(t) = 1 + A+— 7 +\'l^—r 
J + \ > T + M 2 ^ + M 4 



(5.6) 



one may expand the integral in the denominator according 
to 



dydzA(y,z)f+(t) 2 (5.7) 



= a + ai X + + a 2 (X 2 + ± 2A+) + a 3 A + A" + a 4 A^ 2 . 

The numerical values for the coefficients ao-4 as found 
by performing the relevant phase space integrals are given 
in Table |5] We remark that although we neglect isospin 
breaking effects in this subsection, we have employed the 
physical kaon and pion masses [Appendix |A"] . 

Similarly, one can also calculate the dependence of the 
numerator on the form factor parameters A + , A" . In an 
analogous manner to (|5.7(l we write 



P 



'LIPS 



SM 



(5.8) 



60 + &iA+ + b 2 \\ + 63A+ + 6 4 A+A^ + 6 5 A" 2 



iKD 
°0 



bf 



-0.011 ±0.003 



-0.02 ±0.01 



-0.06 ±0.02 



The coefficients bo , 61, and 63 have contributions also from 
the structure-dependent terms, such that we decompose 
them again according to bi = b] B ± 6^ D . b 2 , 64, and 65 
have no structure-dependent part. In our framework, the 
bremsstrahlung amplitude is expressed in terms of a com- 
pletely general (phenomenological) form factor f+(t), 
while the coefficients &f D can be chirally expanded and 
receive their leading contribution at 0(p 4 ). We point out 
that all the coefficients bi depend on the experimental cuts 
E^ ut , O^f, we however suppress this dependence in our no- 
tation. 

We mention that, in principle, the inclusion of struc- 
ture-dependent terms re-introduces a dependence on 
/+(0) by which these terms have to be divided in order to 
arrive at (|5.5(l . However, the uncertainty in the structure- 
dependent terms themselves coming from higher order 
(0(p 6 )) contributions is at least one order of magnitude 
larger than the uncertainty in / + (0), so we do not have 
to worry about a very precise value for the latter. For our 
purposes, we have used the value predicted (parameter- 
free) in one-loop ChPT, / + (0) = 0.977 [22|. 

The values for the coefficients bi can only be found nu- 
merically in this case. Our findings for the "standard cuts" 
E* ut = 30 MeV, 0™* = 20° are collected in Table For 
the values of the low-energy constants see Appendix^ We 
neglect any variation in these constants as we include esti- 
mates of the uncertainties stemming from the 0(p e ) con- 
tributions (see Sect. 14. 2)l that generously cover the range of 
values for Lg, U[ . These uncertainties are quoted as errors 
for the 6f D in Table We describe the precise procedure 
how we estimate these ranges numerically in Sect. Ifi.2l and 
only note for now that the possible corrections are roughly 
30%, which one would naively expect for chiral SU(3). 

The central observation here is that structure-depen- 
dent terms as predicted by ChPT at one loop contribute 
as little as 1% to each of the parameters in Table UH and 
therefore to the total radiative decay rate. 



5.2 Form factor dependence of 1Z 



We are now in the position to give a numerical predic- 
tion for 1Z that depends solely on A + = A + / A^ and A" = 
A"/(A+) 2 . We choose to normalize all parameters by the 
central value A+ = 0.0294, see Appendix^] in order to ex- 
pand in terms of quantities of natural order of magnitude. 



Note that X\_ z is 



a natural scale for A" 



that one would 



obtain e.g. from K* dominance. 
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Table 4. Coefficients for the A+, A+ dependence of TZ , TZ. 
The errors for the Ci are p 6 estimates. All numbers are given 
for the standard cuts. 



n IB (i,o) x io 2 


ci B x 10 3 


c 2 


x 10 4 


c 3 


x 10 4 


0.963 


-0.0 




-1.5 




1.2 


7£(1,0) x 10 2 


ci x 10 3 


C-2 


x 10 4 


C3 


x 10 4 



0.957 ±0.002 0.3 ±0.2 -1.6 ±0.2 1.5 ± 0.2 



the formula 1)5. 5JI and note that the value for Lg used in 
corresponds to A + = 0.032. We have repeated the calcula- 
tion with the matrix element at order p 4 provided in |S] . 
By use of (JEHl and J£3, we find TZ = 0.96 x 10~ 2 , in 
perfect agreement with the value displayed in the fourth 
row in Table 

The important conclusion is that imprecise knowledge 
of the K e z form factor does not preclude a precise predic- 
tion of TZ. 



The numerical prediction is obtained from (|5.5(l , 1)5.7(1 , 
and l|5.8ll . In order to make the form factor dependence 
more transparent, though, we expand TZ according to 

TZ (A+, A'|) = TZ(1, 0){l ± a (A+ - 1) + c 2 (A+ - l) 2 

+ c 3 A'; + ...}, (5.9) 

where we only retain the leading (and most important) 
terms. 

We begin again by considering bremsstrahlung only. 
The numerical results are given in Table again for the 
standard cuts. They demonstrate that TZ lB is extremely 
insensitive to the details of the K e % form factor due to a 
large cancellation of the A + (A" ) dependence in numerator 
and denominator of TZ IB . 4 Furthermore, we note that a 
tree-level calculation of TZ in ChPT would amount to TZ m 
(as there are no structure-dependent terms) with A + = 
A" = (point-like form factors). Numerically, one finds 
T^trec _ o.963 x 10~ 2 , which is even identical to the above 
result in all digits displayed. 

Inclusion of the structure-dependent terms is straight- 
forward from the results given in Table we show the 
numerical results for the complete (IB±SD) coefficients 
also in Table Perpetuating what was done in Tabled 
we again quote errors on all parameters as induced by the 
estimated 0(p 6 ) uncertainties. 

We repeat here the observation made in the previous 
subsection that structure-dependent terms contribute as 
little as 1% to the ratio TZ. In view of the above remark 
about TZ at tree level, the complete one-loop correction 
to TZ is in fact about 1%. Or, to put it even differently, a 
prediction of the radiative decay rate based solely on inner 
bremsstrahlung is expected to have a precision of about 
1%. The parameters Ci are shifted more visibly due to the 
fine cancellation between numerator and denominator of 
TZ, but they remain tiny and do not change the conclusion 
that the form factor dependence of TZ is entirely negligible. 

This is an appropriate place to compare our findings 
with the calculation in 0. We note that there, the branch- 
ing ratio BR(X e 3 7 ) was determined from the chiral am- 
plitude at order p 4 , with E° ut = 30 MeV, (9™* = 20°. As 
is clear from the above, a cancellation of the momentum 
dependence of the form factors does not occur in this case. 
In order to compare with the present calculation, we use 

4 The weaker A+ dependence of TZ was already hinted at in 
a footnote in jl]. 



5.3 Dependence on the experimental cuts 

The near-complete cancellation of all form factor depen- 
dence in TZ suggests the question whether this might be 
accidental due to the specific cuts chosen for the radiative 
decay width. Here, we want to briefly analyze how the 
above findings change when we vary the cuts on E* and 
9* 7 . We restrict ourselves to the bremsstrahlung part of 
the radiative width and TZ. The most important informa- 
tion on the expansion of TZ according to (|5.9I) is collected in 
Table El We find that the coefficients c] B do indeed vary 
considerably for different cuts, but always stay "small" , 
cj B = O(10- 3 ), c™ 3 = O(10- 4 ). The suppression of 

far beyond 10 -3 for the standard cuts turns out to be ac- 
cidental, however. Still, with |A+ - 1| < 0.1, |A^| - 1, 
we find that the K e $ form factor affects TZ at the level of 

io- 4 . 

We should remark here on the latest results for these 
form factor parameters published by the KTeV Collabora- 
tion |2Sj- For the first time, they find significant statistical 
evidence for a non-zero quadratic term in /+(£), together 
with a sizeable reduction of A + . Converted to our conven- 
tions, the combination of their quadratic fits to K e ^ and 
corresponds to A+ = 0.70 ± 0.06, A'j = 1.85 ± 0.40. 5 
While the deviation for the two individual parameters is 
quite sizeable, we find from Table that a simultaneous 
reduction of A + and an enhancement in A" still only mod- 
ify TZ at the order of a few parts times 10~ 4 at best. Our 
finding that TZ is independent of the details of /+ (t) to a 
very large extent therefore remains valid. 



5.4 Isospin breaking 

We have seen above that the ratio TZ can be predicted to an 
amazingly good precision of less than 1%, using ChPT to 
one loop for the structure-dependent terms plus a rough 
estimate of the size of higher-order corrections. At this 
level of precision, isospin breaking corrections - generated 
by real and virtual photons, and by m u — md - become 
relevant, and we discuss these here. 

As soon as one includes virtual photon corrections, one 
also has to take care of additional soft photon radiation in 

5 Similar trends were already noted in the theoretical fits 
in |26|. Note however the latest experimental results from 
NA48 [57], where a free quadratic fit leads to A + = 0.95 ±0.08, 
A" = 0.23±0.52, completely consistent with our central values. 
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Table 5. Coefficients for the A+, A" dependence of 1Z with 
variation of the experimental cuts on EZ , 9^. 





ncut 


ft IB (l,0)xl0 2 


ci B xlO 3 


cP x 10 4 


cP x 10 4 


20MeV 


20° 


1.297 


-2.1 


-0.4 


-2.0 


30MeV 


20° 


0.963 


-0.0 


-1.5 


1.2 


40MeV 


20° 


0.743 


2.1 


-2.6 


4.5 


30MeV 


10° 


1.254 


1.7 


-1.9 


3.3 


30MeV 


20° 


0.963 


-0.0 


-1.5 


1.2 


30MeV 


30° 


0.790 


-1.6 


-1.1 


-0.7 



order to obtain an infrared finite quantity, and we there- 
fore clarify the precise prescription as to what is meant 
by the numerator and the denominator of the ratio R in 
1)2. ip . In accord with the experimental situation |12II13) . 
the denominator denotes the inclusive width for Kl — > 
7T ± e T ^ e ("-7), where (ivy) denotes any number of photons 
of arbitrary energy. The numerator is specified in an analo- 
gous manner: experimental measurements of K e3l require 
detection of at least one hard (E* > E™\ 6»* 7 > 6*™') 
photon in the final state, plus an arbitrary number of ad- 
ditional soft or hard photons. A full calculation of the 
0(a 2 ) contributions in R is beyond the scope of this work. 
Instead, we identify some partial contributions to it and 
give an estimate of the remainder. 

Radiative corrections have been evaluated for K e % in 
28 29 30 . Effects from the quark mass difference m u — m c i 
have been in addition taken into account in [2H| • We note 
the following from that investigation: 

1. One particularly pronounced effect is the electroweak 
correction factor to the Fermi constant, G 2 F — ► SewG 2 ^, 
which contains a large short distance enhancement fac- 
tor [311135] tx log Mz/M p such that Sew - 1 « 2.2 x 
10~ 2 . This factor, however, is universal in the sense 
that it applies identically also to the radiative rate and 
therefore cancels in R. 

2. There are electromagnetic vertex corrections and m u — 
md effects that are collected in a shift in /+(0), which 
can still be factored out as in l|5.3|) . The remaining 
corrections have been incorporated in an expansion of 
the integral / in terms of K e3 form factor parameters 
according to l|5.7[) . The authors of |2U have calculated 
the values for all the parameters ao-4 including correc- 
tions of 0(e 2 p 2 , (m u — md)p 2 )- Their results are col- 
lected in Table 6 In [21], the corrections from real 
photon emission were treated slightly differently from 
what is done here: while there was no upper cut on the 
photon energy, the remaining phase space integration 
of pion and electron momenta was restricted to the 
kinematics compatible with K e % phase space. In order 
to agree with the experimental situation for the case 



6 We are grateful to the authors of |28| for providing us with 
the values for 03 and 0,4 which are not included in the publi- 
cation. 



Table 6. Coefficients for the K e 3 phase space integral, includ- 
ing corrections of 0(a,m u — md)- The numbers for ai to 0,4 
are taken from |28|. For oq, see main text. 



a 


ai 


«2 




a a 


0.09412 


0.3241 


0.4475 


3.080 


6.042 



at hand, we have removed this cut, and have modified 
ciq accordingly, augmenting it by 0.57%. 

3. It remains to estimate isospin breaking effects in the 
numerator of R. A source of potentially large radiative 
corrections are electron mass singularities. Because the 
observed photon is hard and emitted with an angle 
6>* 7 > 20° with respect to the electron, we expect that, 
according to the KLN 32 theorem, these may be ab- 
sorbed into a running electromagnetic coupling con- 
stant. In the present case, the initial state is neutral 
- we therefore evaluate the running coupling at the 
pion mass, rather than the kaon mass, a — > a(l + 

\og(M 2 / m 2 )) . [We stick to corrections of relative 
order a here. Evaluating the coupling instead at the 
kaon mass affects the final result for R by about two 
permille. Up to the number of digits displayed below, 
the final number remains unchanged.] We denote the 
remaining relative corrections by A em . We expect them 
to be small, of the order a/ it ~ 2.3 x 10~ 3 . To be on 
the safe side, we increase this estimate by a factor five 
and take A cm = 0.01. 

4. Finally, we note that part of the running coupling is 
absorbed into /+(0) that contains, in the convention 
of |2HI ) an electron mass singularity as well. We factor- 
ize this piece as before, such that the effect of the mass 
singularity in the numerator is reduced, a — ► a(l + 
jf-logfMl/m^)). As for isospin breaking through 
m u —nid, we expect the effects that cannot be absorbed 
into /+(0) to be tiny, and we neglect them here. 



5.5 Final result for R 

Compared to results quoted in Table pJJ our prediction is 
therefore modified by isospin breaking corrections in the 
following manner: the central value is reduced by 0.2% 
due to corrections in the denominator. We use in the nu- 
merator the running coupling as discussed above, and add 
an uncertainty of ±Z\ Gm . We finally find 

R = (0.96 ±0.01) x 10~ 2 (5.10) 

as our prediction. This may be compared to 1Z IB — 0.963 x 
10~ 2 for bremsstrahlung only, without radiative correc- 
tions. Both R and, as a point of reference, TZ lB are dis- 
played in Fig. pU together with experimental results from 
[TTll 1 2111 311 14| . Note that the corresponding Fig. 4 in [T2"] 
does not properly represent the theoretical result obtained 
in ChPT in [S]: in that reference, the ratio r(K e3l )/ T(K L ) 
was calculated, and an error analysis was not performed. 
In addition, the other two theoretical works OEJ - dis- 
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Fig. 4. Our prediction for R (dubbed "IB+SD") and TZ IB 
("IB") in comparison with experimental results for R from 
NA31 [HI, KTeV 01 E2, KTeV 04 HH, and NA48 [T3]. Note 
that 1Z IB does not contain radiative corrections. All values refer 
to the "standard cuts" E^ ut = 30MeV, 0™* = 20°. 



played with an error bar in that figure - do also not con- 
tain an error estimate. 

We conclude with the observation that the smallness 
of structure-dependent contributions in R precludes a di- 
rect determination of (hadronic) structure effects from R 
alone. In order to extract such effects from experiment, 
one has to resort to differential distributions, which we 
will discuss in Sect.0 



6 Structure-dependent terms: 
Numerical results from ChPT 

In the ratio R, the effect of the structure-dependent terms 
is tiny. On the other hand, in ^2Jj the KTeV Collabora- 
tion attempted to extract two of the structure-dependent 
terms from the E* spectrum. Their result encourages us to 
take up this issue here, in particular so in view of its con- 
nection with the effective theory of the Standard Model. 
The remaining part of this article is devoted to this issue. 

The structure-dependent terms are characterized by 
six amplitudes Vi, Ai. The effect of V3, ^3 is suppressed by 
powers of the electron mass in the decay rate - these am- 
plitudes are only needed for a comparison with the basis 
used in |12j . The Vi , Ai are in general complicated func- 
tions of the variables s, t, and u. At leading order in ChPT, 
they are however real, vary little over physical phase space, 
and may well be approximated by real constants. It thus 
appears that in this approximation, the amplitudes can 
directly be confronted with the KTeV analysis ^21- The 
reason why this is not the case is the following. As we have 
mentioned in Sect. El the IB terms used in the present 
work differ from the ones in |3] [and used by the KTeV 
Collaboration], and therefore the SD terms differ. In ad- 
dition, we use a different set of tensors to decompose the 
SD terms into Lorentz invariant amplitudes. We have dis- 
cussed before why we believe that the decomposition into 



IB and SD terms used here is more appropriate in K e ^ 
than the one originally proposed in . As for the choice of 
a tensor basis, the one used here has the advantage that it 
automatically singles out the two amplitudes whose effect 
in the rate is suppressed by powers of the electron mass. 
The basis used in PITT^] does not have this property, as a 
result of which the interpretation of the various SD terms 
is somehow involved, see below. 

We present the result of our analysis in the following 
manner. First, we discuss numerical results for V^,Aj at 
leading and next-to-leading order in ChPT. We then detail 
the decomposition of the hadronic: tensors V^ U ,A^ V used 
in |12j , and translate our result into the Lorentz invariant 
SD amplitudes used there. 



6.1 Numerical evaluation of the 0(p 4 ) terms 

An important assumption of the analysis of structure- 
dependent terms in ^2] is that these are real and constant 
- which is not really true. Let us therefore investigate in 
what sense real and constant structure functions can be 
taken as reasonable approximations. 

Whereas the leading contributions to the Vj, Ai are 
real, imaginary parts develop at higher orders in the chiral 
expansion. It is shown in Appendix iBl that their effect is 
suppressed in the physical quantities considered in this 
work. More precisely, imaginary parts occur only through 
contributions quadratic in the SD terms and are therefore 
completely negligible here. Concerning the use of constant 
form factors, we have already mentioned that the leading 
contribution to the SD terms indeed is very slowly varying 
over physical phase space. To quantify this statement, we 
average the real part of the ChPT structure functions, 
i.e., integrate over phase space and divide by the phase 
space volume. We quote the standard deviation for this 
average in order to quantify how sensible the assumption 
of the functions being constant is. We use the notation 
(Vi) for the result of this averaging procedure, and quote 
the numbers in units of the kaon mass. 

The numerical results for the structure functions Vi, 
Ai at order p 4 are collected in the first column of Tables 
and 03 at F 2 = F^Fk- [Like for the analysis of R in the 
previous section, we use the central values for Lg, L\ as 
displayed in Appendix lAl and neglect their uncertainties, 
which are generously taken into account through the un- 
certainties that we will attach to higher order terms.] The 
axial terms have only tree-level contributions at this or- 
der and are strictly constant. As for the vector terms, the 
variation in the Vi indicated by the error range in the first 
column of Table is very small. In fact V\ , V<i are domi- 
nated by the counterterm contributions (at a typical scale 
like /i = M p ), which are necessarily constant at 0(p 4 ). For 
comparison, we also show the numerical values for the ap- 
proximations in l|4.3|) in the column dubbed accordingly. 



6.2 Numerical evaluation of the 0(p 6 ) results 

It is desirable to get a handle on the typical size of the cor- 
rections to be expected at 0(p 6 ), and we start the discus- 
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Table 7. Average values for the vector amplitudes Vi in l|3.16[l . 
The symbol (Vi) denotes the average of the real part of the Vi 
in units of Mk- The first column displays the result at order 
p 4 (with the size of variation over phase space indicated) , and 
the second column gives the values in the approximation given 
in 14.311 . The last column contains an estimate of higher order 
contributions, see main text for details. 





o(p a ) 




uncertainty 


(Vi) 


-1.26 ± 0.004 


-1.25 


±0.4 


(Vi) 


0.12 ± 0.002 


0.12 


±0.2 




-0.02 ± 0.001 


-0.02 


±0.1 










±0.1 



Table 8. Average values for the axial amplitudes Ai in 1)3. 8f . 

as given by 0(p 4 ) and C9(p 6 ) ChPT. The central values in 
the second column refer to the order p 6 result at the scale 
H — Alp, with the countcrterms set to zero. For the estimates 
of the uncertainties, see main text. The symbol (Ai) denotes 
the average of the real part of Ai in units of Mk- The term 
A3 was not determined at order p 6 for reasons explained in the 
text. 

Q{p 4 ) Ojp 6 ) 
(Ax) -0.07 ± 0.2 

(A 2 ) -0.30 -0.25 ± 0.1 

(A 3 ) 

(A 4 ) ± 0.4 



sion with the axial terms that we have evaluated analyti- 
cally at order p 6 . Our numerical estimate for these terms 
is obtained by taking their real parts at the scale fi = M p , 
averaged over phase space as before. The contributions 
from the counterterms are the essential uncertainty. We 
have estimated the order of magnitude of these polynomial 
terms in the following manner. The low-energy constants 
depend logarithmically on the renormalization scale |33| . 
It seems unnatural for the constants to be much smaller 
than the change induced by the running of the scale, e.g. 
changing the logarithms by one unit. The shifts in the 
polynomial contributions of A\^ induced by a change of 
the logarithm by one unit is the following, 



Ai >ct 
A 2 , c t 



s + 1 - 2u 
-7t-8u 



}■ 

(6.1) 



In both cases, there are (potentially) large ME- correc- 
tions that could dominate the 0(p 6 ) contributions. As A4 
consists exclusively of a counterterm contribution that is 
scale independent by itself, the above procedure cannot be 
applied here. We use instead an even rougher dimensional 
estimate 

Act = ± „_J!L as ■ (6-2) 



K 



Finally, we do an average of these polynomial terms as 
before and quote the result in the second column of Tabled 
as the final uncertainty at this order. We have not worked 
out the amplitude A3 at order p e because it is only needed 
for a comparison with the amplitudes in |12j at order p 4 
and drops out in the basis used in the present work. 

In order to generate an analogous estimate of the con- 
tributions at order p 6 for the vector terms, a two-loop cal- 
culation is needed. This is beyond the scope of the present 
work, and we content ourselves here with the rough esti- 
mates displayed in the last column of Table These are 
obtained as follows. Concerning Vi, we estimate the con- 
tributions at order p 6 and higher to be of the order of 30% 
of the leading term. As V2 is suppressed at leading order, 



we scale its uncertainty by a factor of 2. Finally, for di- 
mensional reasons, the counterterm contributions to 1/3,4 
are constant. The numbers displayed in the last column 
for are obtained from a dimensional estimate similar 
to the one for A4 discussed above. 



6.3 Predictions for the amplitudes 
used in previous analyses 



In the recent analysis ^2] of radiative K e ^ decays, the IB 
part was taken from It differs from the one displayed 
in H3.13(l through terms of order q and higher, and can be 
obtained from V™ by subtracting these additional terms. 
In addition, a different basis of transverse tensors for the 
SD part was used. In this subsection, we first discuss the 
relation between the Lorentz invariant structure functions 
in the two conventions, and then elaborate on their chiral 
expansion. 

The IB part used in J2| is 



V, 



f = ^ (%>*{/+ ~ 2qWf+} -W v {h- 2qWf 2 }] 

■ (6.3) 



2W„{2pJ+-W v h 



The form factors fi as well as their derivatives f% = dft / dt 
are evaluated with argument t. Here and below, barred 
quantities indicate that the convention from 0, l|6.3(l . 
is used for the inner bremsstrahlung part. The SD part 
changes accordingly, such that the sum of IB and SD re- 
mains the same, 



-SD 



(6.4) 



Four of the eight Lorentz invariant amplitudes were re- 
tained in El,ri21 ancl denoted by A, B, C and D. Here, we 
extend this notation to the remaining four amplitudes and 
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write 



V SD = 



A C 



E 



G 



i c 



K 

^ pCT \Mf pP+ M% P ' P \ qa 

+ * t vp ax p p p'° | (p'q g^~p^q x )^r 
+ (pqgi)-p^ x )i§4 

K 

The relation to the V$ , A4 used here is 
A = M 2 K (y 2 +p'q%) , 



(6.5) 



B = -M 2 K [A 2 + p'WA 3 + M 2 A 4 ) , 
C 
D 



M' K [Vi-V 2 -pqV 3 
M 2 K [a x + A 2 + pWA 3 + pp'A 4 



(6.6) 



M A K (V 3 -V A ) 



E 

F = M%{A 3 -A A 
G 
H 



-m v 3 



-Mk 



where 



Vi = V x 



v 3 = v 3 



2A 2 / + 
p'q 

2A 2 /+ - A 2 / 2 



V 2 = V 2 



-4/+, 
2A 2 / + 



p'qqW p'qqW 
A 2 /, = m- f l (W 2 )-2qWf l . (6.7) 

Equation (|6.6() displays the transformation between the 
basis used in ^2] and in the present work, while (|6.7|l 
presents the changes induced by the difference in the IB 
part. 7 

The above relations allow us to calculate the phase 
space averaged structure functions (A) , (B) , . . . in a 
straightforward manner. To be specific, we use linear form 
factors / + , / 2 , as a result of which only the derivative term 
/+ in V 2 matters. The final result is displayed in Table [5] 
where, for reasons that become clear at the end of this 
subsection, we stick to the values at order p A in the chiral 
expansion. In this approximation, the axial terms are con- 
stant - this is why we do not display an error band in the 
last column in Table We now comment on some basic 

7 In order to check the sign conventions used here and 
in 1^1- 11 21 - where the Pauli metric is used - we have alge- 
braically evaluated the expression of the decay width with 
16.31 1. 16.511 in terms of f+,f2,A,B,C,D and in the limit of 
vanishing electron mass. We found complete agreement with 
the corresponding expressions given in (A1)-(A5) of 0, up 
to an obvious misprint in the line after (A3). The amplitudes 
E—H were not used in |3I12) . 



Table 9. Values of the structure-dependent terms in the KTeV 
conventions, as given by 0(p 4 ) ChPT [with the size of variation 
over phase space indicated]. 



<^> 


-1.34 ± 0.002 


(B) 


0.30 


(C) 


0.08 ± 0.005 


(D) 


-0.30 


(E) 


-0.02 ± 0.001 


(F) 





(G) 


0.02 ± 0.001 


(H) 






features of the choice (|6.3|l for the IB part and l|6.5|l for 
the transverse tensors. As for the impact of the difference 
in the IB part, we note that, expanding in l|6.6l) the form 
factors fi(W 2 ) around q = 0, it is readily seen that the SD 

amplitude V^P indeed differs from the one in the present 
work only by terms of order q and higher, as it must be for 
a reasonable choice of IB. On the other hand, as already 
mentioned in Subsect. 13.21 these additional terms are sin- 
gular at s = M 2 , and can potentially distort the ampli- 
tudes near the boundary of phase space. The difference 
in the choice of the SD part generates more pronounced 
effects. As we have mentioned before, the structure func- 
tions V 3 , A 3 are suppressed by a factor of m^jM 2 , and 
are therefore inaccessible in the electron decay mode. The 
tensor decomposition (|6.5|l does not make use of this fact. 
As a result, certain simultaneous shifts in A, C, E, G (cor- 
responding to a change in V3) or simultaneous shifts in B, 
D, F, H (corresponding to a change in A 3 ) are unobserv- 
able. Measurable combinations are 



M 2 K 
M 2 K 
M 2 K 



M 2 K V 2 , 
-M 2 K (A 2 



MiAi 



D 



H = 



pW 

1V1 K 



E + G 
F + H 



M 2 K {Vx-V 2 ) , 
M 2 K (A 1 +A 2 +pp'Ai) , 



(6.8) 



-M%Vi . 
-M A K Ai 



In other words, the decay width can be expressed in terms 
of the quantities on the left hand side of Q6.8|l . We con- 
clude that e.g. the structure functions C, D - or any lin- 
ear combination thereof - are not measurable in K e3ll as 
long as V 3 ,A 3 are nonzero. In the following section, we 
discuss this point in some more detail. In particular, we 
will provide an interpretation of the quantities C and D 
determined by the KTeV Collaboration in |12| . 

Finally, coming back to Table El we note that, because 
A, B, C, and D are not observables, it does not make much 
sense to work out their numerical magnitude at order p 6 . 
On the other hand, their value at order p 4 will be of use 
in the following section. 
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7 Structure-dependent terms in differential 
rates 

7.1 E* distribution: theory 

Of the various differential rates one may consider, the dis- 
tribution dF I dE* stands out for the purpose of extract- 
ing information on structure-dependent terms, as E* is 
the very variable to distinguish bremsstrahlung and the 
structure-dependent part of the amplitude. In our investi- 
gation, we shall neglect the terms coming from the square 
of the structure-dependent amplitude T SD . Furthermore, 
we make use of the observation made in the previous sec- 
tion that in the one-loop approximation, these structure 
functions are constant to rather high accuracy: we replace 
them in the expression l)B.l(l for the square of the matrix 
element by the averages (Vi),(Ai). We then obtain the 
following decomposition of the photon spectrum: 



dr dr } 



IB 



dE* 



dr v . 



(Ai, 



dE* ^ V dE; " dE^ 
+ o(\T SD \ 2 , AVi, AA, 



(7.1) 



The quantity dEy i /dE; denotes the part of the spectrum 
which is proportional to (Vi), and analogously for 
dT^/dE*. [Remember that we define (Vi), (Ai) to be di- 
mensionless.] The quantities AVi, AAi stand for the errors 
introduced by this approximation. 

In the following, we shall neglect the effect of Vi and 
A4. 8 The objective is to study the distributions dFyJdE;, 
drAi/dE^ in order to quantify the possibility to extract 
(Vi) and (Ai) from data. In order to obtain experimental 
information independent of the measurement of the total 
rate, we follow the strategy of ^2] and only discuss spec- 
tra with arbitrary normalization. Furthermore, we follow 
the procedure in that publication and deviate here from 
the "standard cuts", instead we use 8^ = 5°. We have 
found, though, that such a reduction of the angle cut only 
increases the size (and therefore the expected statistics in 
an experiment) of the bremsstrahlung and hardly has any 
effect on the structure-dependent spectra. The relevant 
photon spectra are shown in Fig. Note that the brems- 
strahlung distribution is scaled down by a factor of 200 
relative to the structure-dependent parts. We observe the 
expected fall-off of dTie/eLE* oc 1/E; as well as the lin- 
ear rise of all structure-dependent spectra for small photon 
energies. As phase space bends them down to zero at maxi- 
mum photon energy, all structure-dependent distributions 
show a maximum (a maximum and a minimum in the case 
of A 2 ), which for V\, V 2 occurs around E* = 80 MeV, for 
Ai slightly higher, around E* ( = 100 MeV. Although the 
Ai spectrum has a form distinct from all others, its mag- 
nitude is far too small to be observable. In view of the 



8 We have verified that the distributions for V4 and A± are 
indeed considerably smaller than the ones discussed here, in 
addition to the fact that both (V4} and (^4) vanish at leading 
chiral order. This holds for all differential rates discussed here 
and in Sect. 01 




200 
E T * [MeV] 

Fig. 5. Photon energy distributions from inner bremsstrahlung 
as well as the various structure-dependent terms. The notation 
d/x / dE^ for the various X refers to 117. lfl . The normalization 
factors are Ny t , A t = 200 iVm = W 3 M K /J\K e3 ). We only cut 

QCUt 



on the electron-photon angle, 6% 



5° [H]. 



chiral 0(p 4 ) prediction Ax — 0, this means that no effects 
of the chiral anomaly are likely to be extracted from the 
photon energy spectrum. 

The remaining three structure-dependent spectra are 
remarkably similar in shape, if not in height. If we assume 
that the experimental accuracy is not sufficient to observe 
the slightly shifted positions of the maxima in the three 
spectra, we have approximately 



dri 



dE* 



2.6 x 



dEy, 
~d~E* 



2.4 



dr Al 

dE* 



(7.2) 



where we have taken the height of the peaks as the mea- 
sure for the proportionality factors, irrespective of the ex- 
act energy where they occur. [In case that more accurate 
data is available, it would be straightforward to incorpo- 
rate a more refined representation of the photon spectrum 
than the one proposed here.] 

Equation (|7.2|l is the main result of our investigation 
of the photon spectrum: 

1. To good approximation, the photon energy spectrum 
originating from the bremsstrahlung amplitude is dis- 
torted by one single function f(E;). The information 
on the SD terms is contained in the effective strength 
(X) that multiplies /(£*), 



dr 
1e 



d^^ 



IB 



dE: 



(X) f(El) , 



7 7 

(X) = (Vi)+0.4(F 2 )+0.4(Ai) 



(7.3) 



The three amplitudes V\, V2, A\ differ mainly in terms 
of the weight with which they contribute to (X) . The 



J. Gasser, B. Kubis, N. Paver, and M. Verbeni: Radiative K e s decays revisited 



15 



latter can be calculated in ChPT, 

' -1.2 0(p 4 ) 



(X) = 



-1.2 ±0.4 0(p 6 



(7.4) 



Note that the uncertainty in the contribution from the 
vector channel has only been roughly estimated here. 
3. In order to measure (X), one may use the represen- 
tation ((7.3(1 for the spectrum, insert the explicit form 
of f(E*) and do a fit to the data with (X) as a free 
parameter. Alternatively, one may take any of the am- 
plitudes Vi.2 or A 1: take it to be constant over phase 
space, and perform a fit to the bremsstrahlung spec- 
trum. The result will be the same. However, it is clear 
that in this manner, one has not determined the cho- 
sen amplitude to perform the fit, but just the effective 
strength (X). 



7.2 E* distribution: experiment 

We now discuss the result of the KTeV analysis ^2 m 
light of the previous subsection. First, we note that in |12| . 
all SD parts were set to zero, except the amplitudes C, D, 
that were taken to be constant over phase space. This 
amounts to the procedure mentioned in point [3] above, 
except that two amplitudes have been retained in |12| . 
while one is sufficient to measure (X). Indeed, finds a 
strongly eccentric error ellipse constraining the parameter 
space for these two structure-dependent terms. In order 
to compare the KTeV result with the above representa- 
tion of (X), we translate the KTeV amplitudes into our 
conventions. We assume a linear form factor /+, use the 
relations (|6.6ll and find that, with A = B = , 



V x = C/M 2 K , V 2 = 4/+(0) 



A l = D/M 2 K 



(7.5) 



The Vi,Ai not listed are zero. In other words, the ampli- 
tudes 1(7.5(1 result in the same photon spectrum as the one 
generated by the amplitudes used in We therefore 
conclude that the effective strength (X) is given in this 
case by 



(X) = C + tan (23°)L> + 1.5M|/+(0) : 



(7.6) 



where we have dropped the bracket notation for C, D, be- 
cause (C) = C for constant amplitudes, and the angle is 
introduced for easy comparison with |12| . The structure of 
1)7.6(1 has been confirmed by the observation made in 1 2 
that it is 

C = cos (25.8°) [C + tan (25.8°)£>] (7.7) 

which is best constrained by the data, with 

C" = -2.5 ti.o(stat) ± 1.5(syst) . (7.8) 

This may be compared with the calculation in the frame- 
work of ChPT. Using l|7.6|l and l|7.4|) . and neglecting the 
small difference in the angle, we find 



C" 



-1.6 ±0.4 



[ChPT] 



(7.9) 



which agrees with 1)7.8(1 rather well. 

While an interpretation of the KTeV result 1(7.8(1 as 
a measurement of the effective coupling is sound, it does 
not allow one to draw conclusions about the size of the 
SD terms themselves because, as we have shown in the 
previous section, C and D are not observable amplitudes 
as long as the amplitude V3 is not negligible. In addition, 
the assumption A = B = made in the analysis of |12j . 
on the basis of the soft kaon approximation, is incorrect 
and invalidates such an interpretation of C even for a 
negligible V3. Chiral perturbation theory may be used to 
illustrate this point: we consider the amplitudes at order 
p A and disregard the structure function V3 altogether, then 
from ((6.6(1 . we find 



(X) 



1.4 (A) + 0.4 {B) + (C) + 0.4 (D) + 1.5 M£ /+(0) 

-1.9 



+0.1 



+0.1 



-0.1 



+0.6 



-1.2 



(7.10) 



where we have again used the phase space average for the 
structure functions, because they are not constant in this 
case. Equation ((7.10|) shows that the main contribution 
to the effective strength (X) is due to the amplitude A, 
while C plays a minor role, and the contribution from D is 
canceled by the one from B. Therefore, the approximation 
of setting A and B to zero is not valid and, consequently, 
(X) is not dominated by (C) + 0.4(1?), and should not be 
taken as a measure of this combination of amplitudes. 



7.3 Systematic errors 

We now discuss one potential source for systematic errors 
in this procedure of determining the effective strength (X) 
and start with the observation that the analysis obviously 
requires a rather precise knowledge of drm/dE*. As we 
consider unnormalized spectra, we are insensitive to over- 
all coupling constants and /+(0), but we should investi- 
gate whether a shift in the K e % form factor parameters 
A+, A" can simulate a contribution to the spectrum sim- 
ilar to the structure-dependent effects. For this purpose, 
we expand a general bremsstrahlung spectrum with arbi- 
trary form factor around our choice for these parameters, 



dn 



IB 



IB 



dEl 



dEl 



A+=l 
A"=0 



(An 



dEl 



dEt 



(7.11) 

The two spectra / dE* and dr^n / dE* are displayed 
in Fig. El with a solid and a dashed line, respectively. We 
observe that they are rather sizeable, but very similar in 
shape to the overall IB spectrum. Only a fine-tuned can- 
cellation of both can lead to a peak-like structure, which 
we find for 



1-Ah 



0.080 ±0.005 



(7.12) 
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N x dr x /dE; 

0.20 



0.15 




0.10 



0.05 



0.00 

200 

e; [Mev] 

Fig. 6. Photon energy distributions from inner brems- 
strahlung, proportional to (A+ — 1) and A" , as well as a fine- 
tuned difference between the two. The normalization factors 
are given by Nax = 100 N- x+ = 10 JV X » = 10 3 M K /r(K e3 ). We 

cut on the electron-photon angle with 0™* = 5° |12| . 



Stated differently, a has to be within this narrow range in 
order for the combination 



dr AX 

dEt 



dE* a dEt 



(7.13) 



to have a maximum. This means that, for example, a si- 
multaneous reduction of A + by 8% as compared to the 
central input value and the introduction of a quadratic 
term in the form factor with A" = (A+) 2 (as suggested by 
K* pole saturation) mimics a structure-dependent contri- 
bution with a peak at roughly similar energies as f(E*). 

The distribution H7.13[l is displayed in Fig. El as a grey 
band at a = 0.08. Note that the strength of the peak is 
not big: for the chosen combination, it is about 10% of the 
dominant spectrum (V\) dry 1 / dE* . To illustrate potential 
effects of this "background" , we compare these findings 
to the latest KTeV form factor measurements j^l], A+ = 

0. 70 ± 0.06, A'| = 1.85 ± 0.40. 9 _Taking into account the 
correlation between A+ and X'l , we find that they lead 
to a = 0.16 ± 0.01, and we conclude the following: 

1. Although these values for A + , A" are very different 
from our assumed central ones, they do not lead to a 
peak-like structure. 

2. Even in the worst possible case with a « 0.08 and 
A" ~ 2, the value for {X) based on the assumptions 
A + = 1, A" = is less negative than the true one. In 
other words, the modulus of (X) would be even bigger 
in the real world, by (20-25)%. 

A more detailed analysis of this background phenomenon 
ought to be performed on real data. 

9 Note the conflicting results in |27) . 



7.4 Other distributions 

We have emphasized that the study of the photon energy 
spectrum, at least with the currently achievable statistics, 
seems to give access to only one specific linear combina- 
tion of structure-dependent terms, which is most sensitive 
to Vi . Ideally one would find alternative distributions that 
are more sensitive to the other terms V2, A\, A2 in order 
to achieve a complete decomposition into the four (main) 
structure functions. The strategy for studying the vari- 
ous possible differential rates is to find those structure- 
dependent contributions that differ in shape from inner 
bremsstrahlung and, in contrast to dr/dE*, from each 
other. We have studied differential rates with respect to 
the other four independent variables E*, i?*, x, W 2 , but 
also to related variables s, t, u, cos0* 7 , cos0* 7 . Where ap- 
plicable, we have used the cuts ££ ut = 25 MeV, 6»™* = 5° 
in analogy to the procedure in |12|. 

One general feature can already be seen from the 
dr/dE* plots and appears in almost all distributions: the 
relative importance of V\ , V2, A\ , Ai is the same in most 
cases, as the integral over a differential rate has to be 
the same, no matter what kinematical variable is studied. 
Therefore distributions tend to be most sensitive to V%, 
followed by V2 and Ai at roughly equal strength. The ex- 
ception to this rule is A2 that shows a sign change in most 
distributions, but again in most cases it is suppressed with 
respect to the other structure-dependent terms by at least 
one order of magnitude. 

We shall only discuss those differential rates in some 
detail that seem to have interesting features. The distri- 
butions in W 2 , cos9* 1 seem to offer no promising pos- 
sibilities to extract information on any of the structure- 
dependent terms as their distributions are too similar to 
the dominant bremsstrahlung one, while those in E*, x, 
or s show A2 curves that have interesting shapes (usually 
with an additional zero), but are probably far too much 
suppressed. 

More promising seem to be the partial rates dr/dE* 
that are displayed in Fig.0 There is no divergent behavior 
visible in these distributions, all of them vanish at mini- 
mum and maximum pion energies, and the partial rates 
for bremsstrahlung as well as for Vi, V2, and A\ have one 
peak in the spectrum (A2 is nearly completely suppressed 
here) . We observe that the bremsstrahlung distribution is 
peaked at high pion energies (for E* — M n w 100 MeV), 
and so are the V\ and A\ partial rates, even though their 
respective peaks occur a bit lower. Distinct from all these 
is, however, the V2 contribution to the pion energy distri- 
bution that is peaked at small pion energies. Although the 
overall sensitivity is again by roughly a factor of 3 smaller 
than that for Vi, this partial rate might be a window to 
access information on the structure function V2. 

We remark, though, that this extraction might again 
be obscured by uncertainties in the form factor / + : the dis- 
tributions dr^ + /dE* and dr^/dE*, defined in complete 

analogy to (jT. 1 1|> . also turn out to be peaked for lower 
pion energies than the total bremsstrahlung distribution. 
Of course this problem can be remedied by more precise 
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N y dr y /dE* 




Fig. 7. Pion energy distributions from inner bremsstrahlung 
as well as the various structure-dependent terms. The nota- 
tion drx/dE^ is chosen in analogy to 1)7,1^ , The normaliza- 
tion factors are N Vi , A t = 200 Mb = W 3 M K /r(K e3 ). The cuts 
= 25 MeV, 6™? = 5° were applied. 



N y dr y /dcose* 
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Fig. 8. Distributions with respect to cos d* £1 from inner brems- 
strahlung as well as the various structure-dependent terms. 
The notation dFxfd cos #* 7 is chosen in analogy to 17. II . The 
normalization factors are N Vi ,Ai = 200 Mb = 10 4 /r(K e3 ). 
The photon energy cut E^ ut = 25 MeV was applied. 



form factor data as provided e.g. in [251127) . A more de- 
tailed study should be performed with actual experimental 
data. 

There is special interest in finding a partial rate with 
a more pronounced contribution from Ai for the following 
reason: as discussed in Subsect. 16.11 this is the only non- 
vanishing contribution of the WZW-anomaly term, while 
A\ vanishes at 0(p A ). We have commented before on the 
sign change in the Ai distributions that often leads to 



cancellations. Fig. |H1 shows a partial rate in which Ai is 
relatively prominent: its contribution becomes relatively 
strong in dT '/ 'd cos 0* in backward direction. The slope 
of the total structure-dependent distribution in backward 
direction, which can be thought of as the second deriva- 
tive with respect to cos 6** 7 at cos #* 7 = — 1 , is potentially 
dominated by A2 . It seems therefore that if effects of the 
chiral anomaly should be visible at all, it might be accessi- 
ble in the distribution with respect to the electron-photon 
angle, in backward direction. 

To conclude this section, we emphasize that this study 
of possible additional partial rates is by no means meant 
to be exhaustive. In particular, certain effects may only be 
visible in double differential rates etc. We defer any such 
more extensive study until experiments give hints about 
the statistical feasibility of these various suggestions. 



8 Conclusions and outlook 

In this paper, we have analyzed various aspects of K e ^ 7 
decays. 

1. In the absence of radiative corrections, the decay am- 
plitude may be decomposed into an inner brems- 
strahlung part (IB) and a structure-dependent part 
(SD). Our construction of the bremsstrahlung ampli- 
tude guarantees that the SD part is regular in the 
Mandelstam plane, aside from the branch points re- 
quired by unitarity. Structure-dependent contributions 
can be parametrized in terms of eight structure func- 
tions Vi , Ai , i = 1, 4. 

2. We evaluate the expression for the width with massless 
spinors. In other words, the electron mass is set to 
zero in the numerator of the relevant terms. In this 
approximation, the contribution from the IB part can 
be written entirely in terms of the K e s form factor /+ . 
Furthermore, the structure functions V3 and A3 cancel 
out. 

3. If this IB-SD separation is applied to the chiral 0(p 4 ) 
representation of the K® 3 ^ decay amplitude provided 
in [S], one obtains leading-order chiral predictions for 
the structure functions (they vanish at order p 2 ). The 
axial terms are constant and given in terms of the 
WZW anomaly, while the vector terms receive con- 
tributions both from loop graphs and the low-energy 
constants Lg and L r w of the chiral Lagrangian at order 
p 4 At this order, all cuts in the loop functions lie 
far outside the physical region, such that also the vec- 
tor terms can be approximated to good accuracy by 
constants. 

4. In order to obtain control of higher order corrections, 
we have analyzed 0(p 6 ) contributions to the structure- 
dependent terms. We have performed a complete cal- 
culation for the axial terms. For the vector ones, we 
have determined the Li x Lj contributions at order p 6 , 
and provided a very rough estimate of the remaining 
diagrams. 

At this order, cuts appear in the physical region, both 
in the axial and in the vector structure functions. The 
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corresponding imaginary parts generate T— odd con- 
tributions in some of the decay distributions. On the 
other hand, in the cases considered in this work, they 
drop out. 

The effect of the cuts on the real parts is diminished 
by the fact that they appear as P-wave rescattering 
(axial) or only in a tiny corner of phase space (vec- 
tor). Dominant uncertainties arise from correc- 
tions, generated by counterterms at 0(p 6 ). 
5. The most precise and stable theoretical prediction can 
be given for the ratio R of the radiative K e z decay 
width relative to the non-radiative one. This ratio 
turns out to be very insensitive to the details of the 
form factor / + , such that the purely hadronic result is 
very precise. Structure-dependent terms yield only a 
1% correction to the bremsstrahlung, such that even a 
sizeable uncertainty in the former affects the precision 
of the total value only at the few permille level. Our 



prediction for E° ut = 30MeV, 6£f = 20°, 



R= (0.96 ±0.01) x 10" 



(8.1) 



deviates from the KTeV results |14U12j , but agrees well 
with the recent measurement of the NA48 Collabora- 
tion [T3j, see Table □ 
6. We have investigated the possibility to measure SD 
terms. We find that the bremsstrahlung spectrum is 
modified by the SD terms essentially by one single 
function /(£"*), and that the different structure func- 
tions contribute with different strength to the effective 
coupling multiplying f(E*). The KTeV analysis [H| 
confirms this observation. In their language, the effec- 
tive coupling is obtained from the combination 

C = cos (25.8°)C + sin (25.8°)L> (8.2) 
of amplitudes C, D, with ^2] 

C" = -2.5 ±i;o(stat) ± 1.5(syst) . 
The calculation in the framework of ChPT gives 



C 



-1.6 ±0.4 



[ChPT] 



(8.4) 



We have shown why the result l|8.3|) should not be 
interpreted as a measurement of the amplitudes C, D, 
but rather as a measurement of the effective coupling 
of the SD terms to f{E*). 
7. We have discussed alternative distributions over phase 
space in a qualitative manner. In order to distinguish 
the vector functions V\, V2, the distribution in pion 
energies might be used. Effects of the chiral anomaly 
are highly suppressed in most distributions. It might at 
best be accessible in the differential rate with respect 
to the electron-photon angle in backward direction. 

Most extensions of this work will depend on the further 
interplay between experimental accuracy and theoretical 
desirability: for example, a complete calculation of the ra- 
diative corrections would be desirable. We have refrained 
here from comparing to the latest KTeV results on R with- 
out cuts on the photon-electron angle 14,, as this also 



would necessitate special care concerning radiative cor- 
rections; this will be considered elsewhere. As indicated 
in the above, a more detailed study of how to disentangle 
the various structure-dependent terms would be possible 
once the extraction of such terms from experiment be- 
comes feasible. 

The most imminent extension of this work, however, 
is to provide predictions also for the other channels. 
An analogous study of the charged channel K^ 3l is most 
straightforward and will be performed in due course. The 
muon channels might in principle lend easier access to 
structure-dependent contributions. The KTeV experimen- 
tal determinations of the ratio R for Kl — > ir^^v^ |14j . 
with improved uncertainty with respect to the previous 
NA48 results [33|, can already be considered as an inter- 
esting starting point for a more comprehensive study of 
radiative decays. 
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(8.3) A Notation 



We denote the charged pion and neutral kaon masses with 
M w and Mk, respectively. In numerical evaluations, we 
use 



M K = 497.67 MeV , M w = 139.57 MeV , 
m e = 0.511 MeV , F n = 92.4 MeV , F K 

The K e 3 form factor is parametrized by 



/+(t) = / + (o) 



1 + a+ ^ + a "m 



(A.l) 
1.22 F„ . 



(A.2) 



As explained in the main text, the precise values of /+(0) 
and A + do not matter in the present context. For nu- 
merical evaluations, we use the parameter-free one-loop 
result El 

/+(0) = 0.977 , (A.3) 



and a central value A^j 
stants we take 



L r 9 (M p ) 



6.3 x 10" 



0.0294. For the low-energy con- 
L r w {M p ) = -4.9 x 10" 3 . (A.4) 
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Lg was chosen such that the chiral one-loop representation 
for A + reproduces A^j_ = 0.0294. The sum Lg + L r l0 is 
then fixed from 7r e 2 7 decays. We express the low-energy 
constants in the following, scale-independent form 0: 



Table B.l. Prefactors that multiply the en , bi etc. 



1 MlM^Ml 

L9 = L ^-5l2^ l0g ^fi ' 

1 M%MtM? 
Ll0 = L ^ )+ 512^ log M b 



(A.5) 



Again, the precise values of Lg and L\g do not matter 



B Traces and decay widths 

Here, we give the explicit expression for the sum over spins 
in |T| 2 in the limit where the relevant traces are evaluated 
at m e — 0, and comment on the relation between K e ^ and 
-fTg3 7 decays in the presence of T-odd terms. 



B.l Traces 

We write 

N- 1 £ |T| 2 = ai f + {tf + a 2 f + {t)5f + + a 3 6fl 

spins 

4 

+ [( b * ReV * + $ RcA l )/+(i) 

8=1 

+ (c.Re^+cfReA^/. 

4 

+ e Yl [( rf < lmV * + d * Im ^)/+(*) 

i=l 

+ (eilmVi + eflmAi)^/ 4 
+ OiV^AlKA,), 

£ - M K 3 q ■ (p' X p e ) , 

su = mUqW)- 1 [/+(*) - U(W 2 )] , 
N = 16TraG% | V US \ 2 M 2 K 



(B.l) 



with 



(B.2) 



With this convention for TV, the right hand side in (jB.lj) is 
dimcnsionlcss. In the limit m e = 0, we immediately have 

&3 = 4 = c 3 = 4 = d 3 = 4 = e 3 = e 5 3 = . (B.3) 
We use the abbreviations 
z pp' — a , zpq — b , z pp e = c , z pp v = d , 
zp'q = e , z p'pe = / , S p'Pi/ = 3 ) zqp e = h, 
zqp v = j , zp e p v = fc , zpT'K = Z , zp'W = m , 



(B.4) 



"1 


4/(e 2 ft) 


61 


4/ (eft) 




4/ft 


0,2 


l/(e 2 h) 


62 


4/ (eft) 




4/(eft) 


«3 


4/e 2 


64 


l/(eft) 


04 


2/h 


Cl 


4/e 


ef 


1 


rfi 


A/ {eh) 


C2 


1/e 


d 


4n/e 


di 


A/ {eh) 


C4 


1/c 


as 


1 




2/ {eh) 


df 


1/ft 


ei 


4/e 


e\ 


1 


dl 


4/(eft) 


e*2 


4/e 


el 


4n/e 


d\ 


2/h 


e 4 


4n/e 


el 


1 



and decompose all the coefficients according to a; = di a^ 
etc., where the prefactors di , &i . . . are collected in Ta- 
ble We obtain the following expressions for the coef- 
ficients at, bi and so on: 



z qW = n, z M 2 = r , 



z 



Mk 2 , 



ax = 2bde (e + /) - e [ft (2 ad - g) + j (e + /) 

+fc(e + 2/)] + ftfcr + 2cd(e 2 + 2e/-ftr), 

a 2 = 4 e {e fc [2 d (2 c - ft) + ft - j (1 + 2 c) - 2 fc 
+2 6 (d - ft + fc)] - 2 e ft Z (d - j) 
+2 ft m [d ft + c (j - 2 d) + fc (1 - b)] } 
-4n {2 6 de / + e ft [2 d (-a + e + /) + g] 
-efj-2ek(f + ah) + 2hkr 
+2c{2def + egh-2dhr)}, 

a 3 = -2 e 2 fc [2 (c - ft) (d - j) + fc (2 b - 1)] 
—2 en |e fc (c + d — ft — j — Z) 
+m [d ft + c (j - 2 d) + fc (1 - 6)] } 
+?i 2 [2e(d/ + C5-afc) + r(fc-2cd)] , 

61 = e [d ft (e + 2 /) + .g ft (6 + c) + j (c / - a ft) 

-fc (0 / + aft)] - ftr (dft + cj - 6fc) , 

62 = j (cefc — e ft Z — eft to) + 6 fc [e (ft — fc) + ft to] 

+n ft (c g — a fc) + d ft [e fc — ft m + n (e + /)] , 
64 = 2e{cefc(e-2.g) + e [-6fc(/ + .g) 

+a fc (j - ft + 2 fc) + Z (g ft + / j - e fc)] 
+CTO {2 g h — ej) + fcrn- (6e — 2 ah) 
+d(e + 2/) (-efc + ftm)} 
+2n{-[e (cg(e-2f)-2bfg 
+a{-ek + 2f{j + k)))] +r[c{-2gh + ej) 
+fc(-6e + 2aft)] + d (e + 2 /) (e / - ftr)} , 

bl = h(de + bg + eg) - j (cf + ah) + k (b f - ah) , 
b\ = c e fc ( j - ft) - 6 fc [ft (e + / - .g) + e fc] 

+ft [a fc (ft — j) + Z (—5 h + ej + fj + ek) — den] , 
bl = fc{-de 2 + a(e- 2/) (j - ft) 
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+b [2 f (g - f) - e (g + /)] } + k 2 (-2 r b + 2 a e) 
+1 [gh(e-2f) + (e + 2f) (/ j - ek) +2rhk] 
+m (deh + bek + 2bfk-2ahk) 
+n [f (de + 2bg- 2a j) - r (bk + dh) + aek] 
+c {k (e 2 - 2 r h + 2 r j) + m (2 g h - 2 f j) 
+r j n + e [2 / k — j m — g (2 k + n)l } , 

ci = e {e k (d — h — j — I) + m (2 hj + b k — dh) 

+c {ek — j m)} — en [g (c + h) + f (d + j) — a k] 
+n r(dh + cj — bk), 

c~2 = 8ek\2hj + bk — dh — cj] 

+4n {ek (c + d - 2 (h + j) - I) 

+m (dh + cj — frfc)} — 4n 2 (d f + eg — a k) , 

C4 = 8e 2 k(df + cg~gh~fj — ak + ek) 

+Aen[e k (/ + g - 2 m) 

+m(-2df -2cg + gh + fj + 2ak)] 

+4n 2 [r (df + cg + k(e-a))~2efg] , 
c\ = 4:[bk(g-f)~ehk + ak(h-j) + ejk 

-ghl + fjl - ghn + fjn] , 
4 = bk(f-g)-k(a-2e)(h-j)+l (gh - f j) , 
cl = 4n[ek(g- f) + r k (h - j) + m (/ j -gh)]. 



where 



(B.5) 

The coefficients for the T-odd terms are 
di = -e f + r h , d 2 = h(f + g + h + j) + e(h- k) , 
d 4 = 2e 2 g + 2eh{f + g)-re(h + j)-2rh(h + j) 1 
4 = 4/ , 4 = h(j-f + g-h)-e(h + k) , 

dl = g(2e + 4f) + 2h{g-f)-r(h+j + 2k), 

(B.6) 

and 

ei = -r n + e (/ + g) , e 2 = 2ek-n(f + g + h + j) , 
e 4 = r (/i + j) - e (/ + .g) > 



5 5 _ 



4(5-/) , 
ef = 4n(/-s) 



e| = / ' - g + h- j 



(B.7) 



B.2 On the relation between Kl and K° decays 

Here we comment on the decay Kl — > 7r ± e^^ e 7 and its 
relation to X — > 7r _ e + i/ e 7 in light of the contributions 
proportional to £ in (|B.1(1 . We neglect CP-violating con- 
tributions and write 



1 



\K L ) - —(\K°)-\K )) 



(B. 



The width for Kl — * 7r ± e =F i'e7 is proportional to 

Jd UPS (Ci+C 2 ) , (B.9) 



d = ]T|T(iY L ^- e +^7)| 2 , 

spins 

C 2 = ^|T(^ L ^^+e-p e7 )| 2 . 



(B.10) 



spins 



In G\ (C2), only the component \K°) (\K )) contributes. 
We use CP to transform the second term to the first one. 
Doing so, all three-momenta of the particles change sign. 
Therefore, terms proportional to £ drop out in the sum 
C\ + C 2 , and the decay width for Kl — > 7r ± e =F ^ e 7 agrees 
with the one for K° —> e + v e it~^f , because in this decay, 
£ drops out as well after integration over the momenta. 
These remarks remain true in the presence of the kine- 
matical cuts considered in the main text, in connection 
with the ratio R. Therefore, up to terms quadratic in the 
structure-dependent terms, only the real parts of the am- 
plitudes Vi,Ai occur in the width and in R. Finally, the 
decay width for K° —> Ti~e + v e coincides with r(K e ^). 
This leads to the expression ()5.5(l for the ratio 1Z. 

For T-odd terms in the context of K* 3 decays, see 



C Invariant amplitudes for K® 3 ^ at order p 



In this appendix, we wish to give a simplified form of the 
if ; °3 7 one-loop amplitudes that is nevertheless as accurate 
as the exact result (that can be found in for all prac- 
tical purposes. As the structure-dependent terms start to 
contribute at 0(q), we intend to retain only terms of or- 
der linear in the photon momentum and neglect every- 
thing that is 0(q 2 ) or higher. In this approximation, all 
the structure functions V1/2/3 can be written in terms of 
the conventional two-point function J(t) plus chiral loga- 
rithms and rational functions of the masses. As remarked 
before, V4 = at this order. We use the following defini- 
tions and conventions: 



Mi = M K , m 1 = M?,. , M 2 = M v , m 2 = M K , 

m; . (C.l) 



Ei = M 2 + m 2 , Ai = M 2 - m 2 , 



the Kallen function 

Ai(t) - \(t,M 2 ,m 2 ) (C.2) 
= t 2 + Mf +mt-2(t (M 2 + m 2 ) + M 2 m 2 ) , 

and the loop functions 

Ji(t) = J Kv (t) . Mt) = J V K(t) , 



Jab(t) = Jab(t) - Jab(0) , 



(C.3) 



Jab(q Z ) = - 



1 r d d i 



ij {2nY(M^ a -P){M 2 ~(l-qY) 
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Our results can be written as follows: 



Vi 



F 2 ' 



(C.4) 



4FH 



2Xi(t) 



tSi - 2M 2 m 2 , m 
— - losr — — 



2,E i + {Ml-Ml) )Ji(t) 
2 t-3Si 



48tt 2 



-O(q) , 



V 2 = — (L g + L 1Q ) 



(C.5) 



1 f M 2 K + M 2 
2F 2 \ t 

J x (t)- ' 



2 - 1 



AM 2 K Ml 



16tt 2 Ai(t) 



16tt 2 



t - (M 2 . + Ml) Ml 
32tt 2 {M 2 k -Ml) ° S Ml 



1 2 

F27 X] 



2M 2 m 2 
-/ -'if^l Ai(t) 

1 v v 7 



(2i + 3 (M| - Ml) + A, 



16tH 



t 



Mfm 2 m; 
16tt 2 A g M? 



3 1 



3(M* - M 2 ) + Aj 
2\i(t) 

x f 2(t - Ei) - A, 



t 2 

Ml -Ml 



2tZ t 
Xi(t) 



10 



{Ml - Ml - A t 



16tt 2 



t-mf + 



M 



K 



M. 



21: 



0(q) , 



^ = ^57 El 3 (( M * - ^ 2 ) t + 6 A 2 ) Ji(t) (C.6) 

i=l 

(327i + (M 2 K - Ml)) (t - S { ) - UMfmj 



2Xi(t) 

X ( J ^-li) 
M 2 m 2 f(M 2 K -Ml)t + 3A 2 



\&TT 2 tA t V 

t-6Si 



Xi(t) 



in- 



327T 2 t 



0(q) 



Note that this expansion necessarily upsets the analytic 
structure as the cuts in the variables t and W 2 coincide in 
the limit of vanishing photon momentum. However, these 
cuts lie far outside the physical region (see discussion in 
Sect. 14. 2|) . Furthermore, despite their appearance, all the 
functions above are regular and smooth at t = and t = 
(Mk-Mk) 2 . The results of the even further simplification 
by setting t = are displayed in the main text (|4.3Jl . 






di 



Fig. D.l. Diagrams that contribute to the anomalous ampli- 
tude Ay,v at order p 6 [in a), the contribution from -K,rj inter- 
mediate states vanishes] . Charges of the mesons running in the 
loops are not indicated. The filled vertices denote a contribu- 
tion from the anomalous Lagrangian at order p 4 . External line 
insertions in the tree diagram of order p 4 are not displayed. 



D Axial form factors at order p 6 



In this appendix, we give the explicit formulae for the 
next-to-leading order corrections to the axial form factors 
Ax, A 2 , and as written out formally in (|4.4[) - (|4.6|) . The 
necessary loop diagrams for this calculation are displayed 
in Fig. ID. II We find the following combinations of loop 
functions and counterterms: 



Sx(s) 
S 2 (s) 
Ti(t) 
T 2 (t) 
Ux(u) 
U 2 {u) 



2H^(s) + ^C r ls s 



16t: 2 



C 2s s , 



H r Kn (t)+H r vK (t) + ^fc r lt t 

T r K At) + T; lK (t) + ^c r 2t t, 

-H r KlT (u)-3H; K (u) + ^C r lu a 



2H^{u) + m; K {u) 



3 

16tt 2 



C r 2u u 



(D.l) 
(D.2) 
(D.3) 
(D.4) 
(D.5) 
(D.6) 



2^tt — Mif ~ Mi) 



16tt 2 



(C^Ml + C{ K Ml) 



X 2 = — 



31 
12 



fJ-7T 



19 
"6 



fJ-K 



+ ^f{C^Ml + Cl K Ml). 



(D.7) 



(D. 
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The loop function H^ b (x) is denned as 



H r ab {x) 



1 



12F 2 \ 

x 



l \(x,Mt,M*) 



Jab(x) + 



24tt 2 



32tt 2 



log 



M 2 M 2 _ x S ab - gjg Mg 
H 4 327rM ab 



-4 



(D.9) 



The other functions can also be written in relatively com- 
pact forms: 



1 



24F 2 



13i 



b( log 



M|M 2 



327r^ \ fi 1 * 
2E K „ + 16A K „ 



M|-M 2 (2Z\^ + t) 



Ml) 



log 



4vr 2 tZ\^ 
(«- 327 X)r )(t- A*, 



log 
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T r nK {t) 
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127T 2 t 
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(D.10) 



32ir* \ yc 1 
2E v k + 8A v k 



9 A 



TjK 
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M 2 ^ 
lo gT7f 
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A K * 



M 2 M 2 K {2A Klv ~t) M 2 
+ 47T 2 tZ\ r , /f ° S M 2 - 

_ (t - 3Er,K)(t - A K «) 
12ir 2 t 
M 2 M 2 
32tt 2 F 2 S /x 2 ' 



JriK{t) 

(D.ll) 
(D.12) 



l/f2 



M 2 , A, 



Af„ 2 



where we have used I7 a j, = M L 

and the two point function J a b(x) as defined in l|C.3|) . The 
combinations of low-energy constants occurring in 14.411 
(|4.6j) and (|D . 1 fl - l|D.8|l are given in Table ID . ll according to 
the numbering in |33j . 



E Inner bremsstrahlung in K® 3 ^ decays 

In this appendix we discuss the separation of the hadronic 
tensor V^ v into an IB and a SD part. To be specific, 
we imagine a calculation in the framework of ChPT to 
all orders and discuss the decomposition there. The rel- 
evant diagrams can be grouped in two classes, displayed 



Table D.l. The coefficients from l|D.lfl - l|D.8fl in terms of the 
renormalized low-energy constants Cj". For example, C{ s = 
4C13 — 10C*JX r + ■ ■ ■• Constants without superscript r are scale 
independent. 
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C23 




C-2u 


C-2-ir 


C-2K 
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24 


-24 
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48 
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-16 
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16 











16 
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-16 


-32 


































48 


































-48 





riWr 
°13 


4 
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-10 
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-8 


-10 





22 


10 


-16 


o 14 


-10 


-10 


-4 


12 


18 


20 


16 





-20 


-16 


48 


^15 


8 


8 


8 


-12 


-12 


-16 


-20 





28 


20 


-32 


/~tWr 
1^19 


2 


2 


2 


-2 


-2 


-4 


-2 





2 


2 





^20 


-2 


-2 


-8 


8 


2 


4 


8 





-20 


-8 


16 


^21 


4 


4 


4 


-4 


-4 


-8 


-4 





4 


4 





^22 


-1 


-1 


-4 


4 


1 


6 


4 


8 


-10 


-4 


8 


n W 
^23 


9 


9 


6 


-6 


-9 


-18 


-12 


-12 


6 


12 


8 



in Fig. IE. II The hatched blobs denote one-particle irre- 
ducible contributions. The diagram b) generates a pole in 
the variable s = (p 1 + q) 2 , at s — M 2 , corresponding to 
the intermediate pion state. We isolate the contribution 
of this pole by writing 



V = 



V^ + ^L [2p u f + (W 2 )-W,f 2 (W 2 )} , (E.l) 
p'q 



where W = p—p' — q. In the following, we assume that this 
is the only singular part at q = in the tensor V^, or, in 

other words, that V^v is regular at q = 0. This is the only 
assumption in the derivation of the final expression for the 
IB term. We have checked that it is true at one-loop order 
in ChPT, see below, and we see no reason why it should 
not be correct to any order, and thus true in QCD. Next, 
we write this regular part as 



%u = v g^ v + vi p'^qy + v 2 W^q u + v 3 p'^Wy 



(E.2) 



K° 




Vf" 



a) 




b) 



Fig. E.l. Diagrams for V^ v , evaluated in the framework 
of ChPT. The hatched blobs denote one-particle irreducible 
graphs. 



J. Gasser, B. Kubis, N. Paver, and M. Verbeni: Radiative K e s decays revisited 



23 



The Ward identity (^.fijl generates three conditions on Vuv, 

v + vi p'q + v 2 qW = 2A/+ - / 2 , 

v 3 p'q + v 6 qW = 2A/+ - A/ 2 , (E.3) 
Vi p'q + V 5 qW = 2A/+ , 

with 

A/i = fi(t) - fi(W 2 ) . (E.4) 

The first equation can be solved for vq. Furthermore, we 
set 

2A/+ 2A/+ - A/ 2 , , 

qW qW 

and are left with 

v 4 p'q + v 5 qW = , 

(E.6) 

i>3 J? g + v e qW — . 

At this stage, we use the fact that the Lorentz invariant 
amplitudes Vi are defined for any value of the kinematic 
variables p'q, qW , and that the amplitudes are assumed 
to be non-singular at p'q = 0. It then follows that are 
proportional to p'q, 

V5,6 = -P'q Vi,3 ■ (E.7) 

where the sign and the numbering is chosen for conve- 
nience. Finally, we obtain 

v 3 A = 9^*3,4 ■ (E.8) 

Collecting the results, we find that can be written in 
the form displayed in (|3.16|l , with 

(V 1 ,V 2 ,V 3 ,V 4 ) = {v 1>V2 ,v 3 ,V4) . (E.9) 

Appendix [U] contains the explicit expression of the form 
factors Vi in the limit q — 0, illustrating the fact that they 
indeed are non-singular at q = at next to leading order 
in ChPT, as mentioned above. 
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